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Abstract 

A prototype experiment to measure the neutron electric dipole moment (nEDM) by 
spin-rotation in a non-centrosymmetric crystal in Laue geometry was carried out in 
order to investigate the statistical sensitivity and systematic effects of the method. 
The statistical sensitivity to the nEDM was about 6 • e-cm per day and can 
be improved by one order of magnitude for the full scale setup. Systematics was 
limited by the homogeneity of the magnetic field in the crystal region and by a new 
kind of spin rotation effect. We attribute this effect to a difference of the two Bloch 
waves amplitudes in the crystal, which is caused by the presence of a small crystal 
deformation due to a temperature gradient. In a revised scheme of the experiment, 
this effect could be exploited for a purposeful manipulation of the Bloch waves. 
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1 Introduction 


The search for a hnite nEDM is a prominent example of the quest for new 
sources of CP violation and thus for physics beyond the Standard Model. The 
most precise experiments today were carried out using the Ramsey resonance 
method and ultra-cold neutrons (UCNs) [1,2]. Further progress is presently 
limited by systematics [3] and the low density of UCNs available. 

In general, the statistical sensitivity of an experiment to measure the nEDM 
is determined by the product Et\/N, where E is the value of the electric held, 
r the duration of the neutron interaction with the held and N the number of 
the counted neutrons. New projects to measure the nEDM with UCNs aim 
to increase the UCN density and thus N by orders of magnitude (see [4] for 
a recent overview). In contrast, experiments with crystals exploit the electric 
held inside matter, which can, for some crystals, be orders of magnitude above 
the electric held achievable in vacuum. 

EDM experiments with absorbing crystals were pioneered by Shull and Nathans 
[5]. The hrst one based on the interference of electromagnetic amplitude with 
the imaginary part of the nuclear one [6]. The hrst who paid attention to 
the presence of a spin dependent term due to interference of nuclear and 
spin-orbit parts of scattering amplitude in the neutron interaction with a non- 
centrosymmetric non-absorptive crystal were Abov with his colleagues [7]. 
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Spin-rotation in non-centrosymmetric crystals due to such interference as a 
way to search for a nEDM was hrst discussed by Forte [8]. The corresponding 
spin-rotation effect due to spin-orbit interaction was experimentally tested by 
Forte and Zeyen [9]. A similar to [8], but more detailed theory of neutron 
optical activity and dichroism for diffraction in non-centrosymmetric crystals 
has been developed in [10]. Authors of the works [11,12] have shown and ex¬ 
perimentally proved that the interference of nuclear and electromagnetic parts 
of the scattering amplitude leads to an existence of a constant strong electric 
field, acting on a neutron during all time of its movement in the noncen- 
trosymmetric crystal. This held was measured hrst [12] in the neutron Lane 
dihraction experiment, the measured value being coincided with the calculated 
one. 

Recently a new method of a neutron FDM search was proposed [13,14,15,16]. 
The value of the electric held in this method was determined experimentally 
for quartz to F ^ 2 ■ 10® V/cm [12,17]. The interaction time of the neutron 
with the electric held is shorter than in UCN experiments and can reach 
T ^ 1 — 2 ms [18,15,19]. The statistical sensitivity of the method prohts from 
the higher hux of the used cold neutrons, compared to UCNs available today. 
In a test experiment [19] we have determined the statistical sensitivity of the 
method and have found that with existing quartz crystals and the hux of the 
PFIB beam line [20] of the ILL one can reach ~ 6 ■ 10“^®e-cm per day, about 
the value of the most sensitive published UCN experiments [1,2]. 

The aim of the experiment presented in this paper was to conhrm the statis¬ 
tical sensitivity with a prototype set-up and to investigate systematic errors 
of the method. 
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2 Laue diffraction method to search for a nEDM 


In non-centrosymmetric crystals the diffracted neutrons are moving in two 
Bloch states exposed to opposite electric helds [12,13] because of the shift of 
the “electric” planes relative to the “nuclear” ones ^. In the reference system 
moving together with the neutron the electric held is seen as a magnetic one. 
The interaction (Schwinger interaction) of the neutron spin with this effective 
magnetic held = \[E^ x v] results in a spin precession around iT®. The 
spin rotation angle for the two states is [13] 




r = ±pn 


eEgL 

9 ’ 

nipC^ 


( 1 ) 


where r = L/v\\ is the interaction time (n|| being the component of the neutron 
velocity parallel to the dihracting planes and L the thickness of the crystal) 
and fi = /ineh/2mpC^ the neutron magnetic moment with pn = —1.9. The signs 
± refer to the Bloch states 1 and 2, respectively. For a longitudinally polarized 
beam with the incident polarization Pq, the hnal neutron polarization P is 
longitudinal and given by [15,16] ^ : 

P = Po cos A4 2 = Po cos ( . (2) 

V TTlpP ) 


It can be decreased down to zero by choosing the crystal thickness Lq such 
that A (;/)®2 = ±7r/2. The calculation for the (llO)-plane of a-quartz gives 

^ In our notation, “nuclear” or “electric” planes are determined by the positions of 
the maxima of the corresponding periodic potential. 

^ This result is obtained by averaging the Pendellosung oscillations over the Bragg 
angles. The angular period of these oscillations in the experiment was 10 ^ rad 
and the angular divergence of the neutron beam 10 ^ rad. 
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Lq = 3.5 cm. 


We consider Bragg angles 9b close to n/2, i.e. v || Pq is almost parallel to 
Eg. For each of the two Bloch states, a component of the polarization vector 
perpendicular to Eg builds up in the XY-plane due to the Schwinger precession 
Eqs. (1,2). These components as well as the electric field Eg have opposite signs 
for the two Bloch states. The interaction of a hnite nEDM with the electric 
held Eg results in a precession of the polarization vector around Eg, thus 
creating a component in Z direction, with the same sign for the two 

Bloch states. jg gjygj^ i^y [ 15 ^ 16 ]; 


pEDM _ dPPgLp _ 4P ^ c ^ 1 

nhv\\ TT/i ncos6*B n/2 —9 b 


(3) 


Here D is the nEDM. Rotating the crystal by the angle 29b ^r) changes 

the sign of Eg and thus that of P^^m jg experimental signature of a 
nEDM. 


The principal scheme of the method is shown in Fig. 1. We assume that the 
environmental magnetic held is low enough as to avoid any further spin rota¬ 
tion. For the two crystal positions right (R) and left (L) with the same Bragg 
angle but with opposite directions of the electric held, the polarizations 
have opposite signs whereas a residual polarization has the same sign. In the 
coordinates of Fig. 1, is given by the diherence of the Z components of 

the hnal polarization vectors for the two crystal positions R and L. High preci¬ 
sion of the crystal rotation and “zero” magnetic held conditions are necessary 
to exclude systematic errors and to select the EDM ehect. 

/^From Eq. (2) follows that the ehect due to the Schwinger interaction does 
not depend on neutron properties such as the energy, the wavelength or the 
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Fig. 1. Principal scheme of the experiment for a nEDM search by the Lane diffraction 
method. The presence of a nEDM leads to a small Z-component of the polarization, 
which has opposite signs for the two crystal positions R and L. 

Bragg angle. It is determined by the property of the crystal and by the fnnda- 
mental constants only. For a given crystal it is the same for all Bragg angles. 
In contrast, the EDM effect Eq. (3) depends on the Bragg angle and grows 
strongly for O-q 7r/2. Thus, carrying out the measurement for two different 
Bragg angles gives an additional way to eliminate false effects related to the 
Schwinger interaction. 


3 Experimental set-up 

In order to avoid systematic effects, the two crystal positions have to be fully 
identical. Any gradient of the residual magnetic held or a temperature gradient 
over the crystal violates this requirement and, hnally, can result in a systematic 
offset on the hnal polarization which could mimic the searched effect. 

The experimental setup is shown in Fig. 2. It basically consists of the part 
responsible for the preparation of the beam with the desired neutron wave¬ 
length and polarization (neutron velocity selector, polarizer, hipper), the part 
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responsible for the zero-field condition on the sample and spherical polarime- 
try (Cryopad with nntators [21]), and the analyzer of the hnal polarization. 



Fig. 2. General layout of the experiment. 

The (110) plane {d = 2.456A) of a quartz crystal was used in the experiment. 
The size of the crystal was 140 x 140 x 35 mm^ and its mosaic ~ 1 angular 
second for the whole body of the crystal. 

The experiment was carried out at the end position of the neutron guide H53 
(instrument PFl) of the ILL. The guide is connected to the horizontal cold 
source and provides a capture flux of 3.5 ■ lO^n/cm^s. The velocity selector 
served to preselect neutrons with a wavelength of (5.0 ± 0.3) A and thus re¬ 
duced the background. To obtain a homogeneous and wavelength-independent 
incident polarization, two super mirror polarizers were used in crossed geom¬ 
etry [22]. With this set-up, an incident polarization of about 99.5% is ex¬ 
pected [22]. The experimental values AP = 97.8% for the direct beam and 
AP = (97.5 ± 0.2)% for the diffracted beams were limited by the analyzing 
power A of the single super mirror analyzer. The neutron spin could be flipped 
by a resonance flipper. 

For spherical polarimetry, the so-called Cryopad-II was used. It consists of 
two concentric cylindrical Meissner magnetic screens which prevent the pen- 
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Fig. 3. Angular dependence of selected elements of the polarization matrix Pij. 

etration of the external field inside the cavity at the height of the beam, fi- 
metal screens located above and below Cryopad-II, together with the Meissner 
screens, assure a reduced magnetic held in the sample chamber. Nutators and 
superconducting coils are used to select the direction of the incident polar¬ 
ization vector and the component of the outgoing polarization vector that is 
being analyzed. The residual magnetic held inside the cavity was of the order 
of a few mG. The cross-section of the beam at the sample was 0.8 x 1.7 cm^. 


The cylindrical geometry of Cryopad leads to a discrepancy of the coordinate 
systems between the crystal positions L and R (see Fig. 2) and limits the 
accuracy of the measurement of the polarization vector at about ±10 mrad. 
We estimate that this factor and the residual magnetic held result in a total 
directional uncertainty of the polarimetry of about 20 mrad. This value agrees 
qualitatively with the typical precision obtained with this second-generation 
Cryopad of 35 mrad [21]. 

Although the superconducting screens are obviously thermically isolated, the 
air temperature in the center of Cryopad was about 0.5 K higher than that 
close to the Cryopad walls. 




4 Experimental results 


The intensity of the diffracted beam for the experimental geometry was about 
3-6 n/s for the “grey” polarization direction. It coincides with the result ob¬ 
tained in our previous experiment [19]. The relatively small intensity was 
mainly limited by the small size of the incident nutator restricting the beam 
cross-section to 0.8 x 1.7cm^. The corresponding statistical sensitivity for the 
nEDM is 6 ■ 10“^^ e-cm per day. This value can be improved by about one 
order of magnitude by using a more intense neutron beam (for instance the 
instrument PFIB [20] of the ILL) and increasing the sizes of the beam and of 
the quartz crystal. 


Examples of the measured angular dependences of the elements of the polar¬ 
ization matrix are shown in Fig. 3. Here, i indicates the direction of the 
incident polarization and j the analyzed component of the final polarization. 
The presence of a nEDM would lead to a difference of Py^ for the two crystal 
positions. The observed difference (triangles in the Fig. 4) is consistent with 
zero and can be used to estimate the nEDM D = (3.5 ± 1.6) ■ 10“^^ e-cm 
(only statistical error stated). This value, obtained within a few hours of data 
collection, is 20 times more precise than the result of the crystal diffraction 
experiment of Shull and Nathans [5]. We did not try to improve the statistical 
accuracy further since we found large hnal polarizations P^y, Pyy, see Figs. 3 
and 4, in contradiction to the expectations. 
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Fig. 4. Difference of selected elements of the polarization matrix Pij between the 
two crystal positions R and L. 

5 Discussion 


To our understanding, the only reason for measuring large elements Pxy, Pyy of 
the polarization matrix is a difference between the two Bloch waves amplitudes 
excited in the crystal. ^ 

A difference of this type can be explained by the Borman effect (different 
absorptions for the two Bloch waves). However, our estimations for the (110) 
plane of a quartz crystal show that this mechanism can account for not more 
than 10% of the observed effect. 

A deformation of the crystal could be another possible reason. It is well known 

^ Note that a nonzero value of the matrix element Pxy indicates a neutron spin 
rotation due to spin-orbit interaction in Laue diffraction. Earlier, such an effect was 
predicted only for absorbing crystals [10]. 
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that the interaction of neutrons or X-rays with elastically deformed crystals 
strongly differs from the undeformed case [23,24,11]. In our experiment, such a 
deformation could arise from the temperature gradient of 0.5 K in the Cryopad. 
The thermal expansion of quartz is AL/L ^ 10“^AT/K. This has to be 
compared with the Bragg width AAb/A ~ 10“^ of the (110) plane of the 
crystal. 


The trajectories of neutrons in a deformed crystal can be described by the 
so called “Kato forces” [25] that are determined by the value of the crystal 
deformation. For a constant gradient in the interplanar distance the neutron 
trajectories inside the crystal are given, see [11,25]: 


d'^z tan^ 6b 

^ = ±- ^ 9 ^: 

oy^ mo 


(4) 


where mo = 2dFg/Vc is the so called “Kato mass” (Fg is the neutron structure 
factor of the reflection, Vc the volume of the unit cell, and d the interplanar 
distance), g = 2Ti/d the reciprocal lattice vector, and ^ describes the crystal 
deformation (d = do(l -|- C^))- The signs ± in Eq. (4) correspond to the two 
Bloch waves. Eq. (4) holds for crystal deformations small compared to the 
Bragg width of the reflection. 

The right part of Eq. (4) is proportional to the square of tands- In our ex¬ 
periment this factor was rather high: tan^ 9b ~ 100 ... 1000 {9b ~ (84 — 87)°). 
This results in a very high sensitivity of the neutron trajectories to even small 
deformations of the crystal. 

We calculated neutron trajectories in a deformed crystal of 14 x 3.5 cm^ for 
the temperature gradient AT = 10“^K/cm and 9b = 86° {ia.v? 9 b = 200). 
The results are shown in Fig. 5. The left plot corresponds to a constant tem- 
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Fig. 5. Examples of neutron trajectories in the deformed quartz crystal ((110) plane, 
d = 2 . 45 A, = 86°, crystal size 35 x 140mm^). The neutrons enter the crystal from 

the left side in the center. Left: Constant temperature gradient d = do{l+^z). Right: 
The center of the crystal has a higher temperature than its sides, d = do(l + ^kl)- 
This was the case in the experiment. 

perature gradient from one side of the crystal to the other and the right one 
to a gradient from the center to both sides. The two cases behave differently: 
For the constant gradient the crystal deformations do not result in different 
amplitudes for the two Bloch waves at the exit surface. For the gradient from 
the center to both sides, one Bloch wave ('h(l) in Fig. 5) is focused, the other 
one defocused. Consequently, the amplitudes of the two Bloch waves at the 
exit surface of the crystal are not equal. The outgoing neutron wave is the 
superposition of the two Bloch waves where the spin was rotated in opposite 
directions. Because of their different amplitudes, the depolarization effect de¬ 
scribed by Eqs. (1,2) is incomplete. For equal amplitudes, the X components 
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of the polarization vectors of the two Bloch waves cancel everywhere (com¬ 
pare Fig. 1). For different amplitudes this is not the case and thus the total 
polarization vector is rotated in the XY-plane. In our model, it is the new spin 
rotation effect which is responsible for the high values of the elements Pxy and 
Pyy observed in the experiment. 

It is important to point out that, although is measured via the Z- 

component of the hnal polarization vector, the surviving polarization in the 
XY-plane can cause an offset to p^^M^ example due to a residual magnetic 
held that turns the polarization toward the Z direction. 

On the other hand, the effect permits to manipulate the amplitudes of the 
Bloch waves by applying a temperature gradient. For the inverse sign of the 
gradient, for example, 'h(l) is defocused and 'h(2) focused. Manipulating the 
sign of the temperature gradient is equivalent to changing the sign of the 
effective electric held seen by the forward dihracted neutrons. By a high tem¬ 
perature gradient one of the Bloch waves can even be fully suppressed. 

To test our explanation of the observed ehect, we intentionally introduced a 
temperature gradient along the crystal by installing a small electric heater 
close to one surface of the crystal and measured the spin rotation ehect {Pyx) 
for diherent heating powers. The results are shown in Fig. 6: Heating the 
crystal with an electric power of about 2W changed the sign of Pyx- The 
process had a relaxation time of the order of half an hour. Reducing the 
power of the heater to 1W resulted in the compensation of the spin-rotation 
ehect for the used crystal position. These results are in agreement with our 
explanation. 
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Fig. 6. Dependence of Pyx on the heating power applied to one side of the crystal, 
for the crystal position R. 

6 Conclusions 

A prototype experiment to measure the nEDM by Laue diffraction in non- 
centrosymmetric crystals was carried out. This test allowed us to determine 
experimentally the statistical sensitivity of the method and to investigate pos¬ 
sible systematic errors. 

The experiment conhrmed our expectation as regards the high statistical 
sensitivity: already the hrst run of a few hours resulted in the value D = 
(3.5 ± 1.6) ■ 10“^^e-cm which is 20 times more precise than the result of the 
previous experiment based on crystal diffraction [5]. With a dedicated instal¬ 
lation, a statistical precision of the order of 6 ■ 10“^® e-cm per day can be 
achieved, comparable to the present nEDM experiments with UCNs and the 
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Ramsey resonance method [1,2]. 


Beside the statistical precision the most important issue is the question about 
systematic errors. The main source of systematic errors is a residual magnetic 
held in the zero-held cavity where the crystal is installed. In our installation we 
used an old version of Cryopad as cavity, which has a relatively high residual 
magnetic held of the order of a few mG. This held limits the EDM experiment 
to a precision of about 10“^^ e-cm. We believe that it is possible to build a ded¬ 
icated spherical polarimeter based on the Cryopad ideas to reach a precision 
of the EDM experiment of the order of a few times 10“^® e-cm. 


The performed experiment also allowed us to discover a new ehect: for for¬ 
ward dihraction of polarized neutrons by a non-centrosymmetric crystal of 
quartz (110 plane) a large (~ 7r/4) spin-rotation ehect induced by a temper¬ 
ature gradient was observed. We identihed this ehect as a combination of the 
Schwinger interaction of the neutron spin with the crystal interplanar electric 
held and the modihcation of the Bloch waves amplitudes by a temperature 
induced deformation of the crystal. 


The consequences of this new spin-rotation ehect are twofold: on the one hand 
it imposes a limitation on the proposed method to search for a nEDM and 
requires serious revision of the method, on the other hand it can serve as a new 
tool to manipulate the sign and the value of the crystal interplanar electric 
held and, hence, the sign of the nEDM ehect. These possibilities open a new 
road to develop a nEDM experiment exploring the high interplanar held of 
non-centrosymmetric crystals. 
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